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\ Abstract: This paper is concerned with optimal time-decay estimates of solutions of 

f-H \ the Cauchy problem to a model system of the radiating gas in W 1 . Compared to Liu and 

Kawashima (2011) |22] and Wang and Wang (2009) [33], without smallness assumption of 
initial perturbation in L 1 -norm, we study large time behavior of small amplitude classical 
solutions to the Cauchy problem. The optimal H N -noTm time-decay rates of the solutions 
in R n with 1 < n < 4 are obtained by applying the Fourier splitting method introduced in 
^-j- \ Schonbek (1980) |27j with a slight modification and an energy method. Furthermore, basing 

on a refined pure energy method introduced in Guo and Wang [6] (2011), we give optimal 
L P -L 2 (R 3 ) decay estimates of the derivatives of solutions when initial perturbation is bounded 



in L p -norm with some p € (1,2]. 

Key Words: the radiating gas; Fourier splitting method; negative Sobolev space; opti- 



^ ■ mal time-decay rate. 



1. Introduction 

In the present paper, we consider the large time behavior of the Cauchy problem for a 
model system of radiating gas, taking the form of 

n 

u t + E fj(u) Xj + divq = 0, x £ W 1 , t > 0, 

(1.1) 

-Vdivq + q + \7u = 0, x G W n , t > 0, 

with initial data 

u(x,0) = u (x), xef. (1.2) 

Here unknown functions u = u(x, t) : W 1 x [0, oo) — > R and q = q(x, t) : R n x [0, oo) — > R n 
represent the velocity and radiating heat flux of the gas, respectively. The notations V and 
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div are the ra-dimensional gradient and divergence. f(u) = (fi(u), ■ ■ ■ , f n (u)) £ K n is a given 
smooth function of u satisfying fj(u) = 0(\u\ 2 ) (j = 1, • • • , n) for u — > 0. 

System (1.1) simplifies the model for the motion of radiating gas in several space vari- 
ables. Indeed, in a certain physical situation, (1.1) is well approximated to the fundamental 
system describing the motion of radiating gas: 

p t + div(pu) = 0, 

(pu)t + div(pu <8> u + pi) = 

(1.3) 

W + ¥ )]t + div[/w(e + l -^f)+pu + q} = 

- Vdivq + a x q + a 2 V6 4 = 0, 

where p,u,p,e,0 are respectively the mass density, velocity, pressure, internal energy and 
absolute temperature of the gas, while q is the radiative heat flux. a\ and a 2 are given positive 
constants depending on the gas itself. The simplified model (1.1) was first investigated by 
Hamer [7], and the reduction of the full system (1.3) to (1.1) was given in [5] \7\ [15], 

A lot of important works have been done on system (1.1). For one-dimensional case, we 
refer to [3] I3CH [3"T] for L 1 stability results, [14 |, I16|, IT9j for a singular limit and relaxation limit, 
[TTlII2lII3liniIIEll2Ill2Il|231|21l|25]for shock waves, El [TO] for diffusion waves and [E] 
for rarefaction waves. 

However, there are fewer studies for system (1.1) in the case of multi-dimensional space. 
Recently, Francesco in [2] obtained the global existence of weak entropy solutions of system 
(1.1) and the relaxation limits. Later, Wang and Wang in [33] investigated pointwise esti- 
mates of solutions to (1.1)-(1.2) by means of a detailed analysis for Green's function. More 
recently, the asymptotic decay rates toward the planar rarefaction waves based on the L 2 - 
energy method are obtained in [3] for 2-dimensions and in [5] for n-dimensions (n = 3,4,5), 
respectively. The asymptotic behavior of solutions to the diffusion waves was studied in 
)22[ 126] , For the related study of decay rates to the problem (1.1)-(1.2), we mention |22| 
and [I] only. [22] studied large time behavior of solutions to the problem (1.1)-(1.2) with 
small initial data in L 1 -norm perturbation by using a time- weighted energy method. In [I], 
under the assumption that ||«o||i,i is bounded, Duan et al. showed that the optimal L 2 -norm 
time-decay rate of solutions is (1 + t)~* , whereas time-decay estimates of the derivatives of 
the solutions have not been considered. 

The purpose of this paper is to establish the optimal time decay rates of the solutions 
and the derivatives of the solutions to the Cauchy problem (1.1)-(1.2) without smallness 
assumptions on the initial data in L 1 -norm. It seems that the usual energy method basing 
on the linearization analysis doesn't work. By the way, our first decay result is inspired 
from Schonbek [281 129] . where the well-known Fourier splitting method is established to get 
optimal decay rate of solutions for the incompressible Navier-Stokes equations in L 2 -norm 
or ilP-norm. In the present paper, we generalize the Fourier splitting method with a slight 
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modification to deal with the problem (1.1)-(1.2) for R n with 1 < n < 4. Our second decay 
result is illuminated by a recent work of Guo and Wang [6], where they developed a new 
method to establish optimal time decay rates of solutions to the Cauchy problem for the 
compressible Navier-Stokes equations and the Boltzmann equation. The main idea was to 
combine scaled energy estimates with an interpolation between negative and positive Sobolev 
norms to get the time decay rate for dissipative equations. By employing this new method, 
we obtain the optimal L P -L 2 (R 3 ) decay estimates of solutions to the problem (1.1)-(1.2). 

First, we cite a temporal global existence result established in [22"1 [33] . 

Proposition 1.1. ([22], Theorem 2.1; [33], Theorem 1.1) Assume that u (x) G H N (R n ) (n > 
1) for an integer N > [8] + 2. There exists a small positive constant 5o such that if 
Eo = \\ u o\\h n (R n ) — ^0) then the problem (1.1)-(1.2) has a unique global solution (u,q)(x,t) 
satisfying 

u G C([0, oo); H N (R n )), Vu G L 2 ([0, oo), H 1 *" 1 ^ 71 )), 
q G C([0, oo); H N (R n )) n L 2 ([0, oo),H N+1 (R n )). 

Moreover, the solution verifies the following uniform energy estimate 

+ J (l|Vn(T)||^_ 1(Mn) + ||g(T)||^ +1(Rn) )dT < CEl (1.4) 

Our main results in this paper can be stated as follows: 

Theorem 1.1. Let 1 < n < 4. Suppose || < 8q -C 1. If further, uq G L (R n ) 

(which need not be small), we have 

\\D l u(t)\\ L2{Rn) <C(l + t)- r i-K 1 = 0,1,.-. ,N; (1.5) 
||£>^)||l2( K «) <c(i + i)- f "^, i = 1,2,--- ,iV-l. (1.6) 
Remark 1.1. One can rewrite the system (1.1) as a decouple system of (u,q)(x,t): 

n 

u t + y~] fj(u) x . = -u + (I - A)~ v 

U (1-7) 

g = -(J — A) _1 Vu. 

To obtain the time-decay rates of q(x,t), from (1.7)2, it suffices to prove the estimates on 
u(x,t), i.e. (1.5). When Z = 0, (1.5) can be obtained directly by using the usual Fourier 
splitting method. In addition, for any integer n > 1, (1.5) also holds true, see Proposition 
3.1 in Section 3. This is consistent with the result in pp. When Z > 1, to prove (1.5), the 
estimates in the Proposition 1.1 and the smallness of Eq should be employed, see details in 
the proof of Proposition 3.2. 
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Theorem 1.2. Suppose ||^o||ijW(R3) < <5o <C 1. If further, uq G H S (]R 3 ) (which need not be 
small), for some s G [0, |), there exists a positive constant Co such that 

IKi)ll^ (R 3) < Co, lk(*)ll^ (R3) < Co, (1.8) 
and the following decay estimates hold: 

\\D l u(t)\\ HN - lm <C (l + t)-^, Z = 0,---,iV, (1.9) 

P'ff(t)ll^-i-«(ms) ^Coa + t)-^, / = (),••• ,JV-1. (1.10) 

By employing the Hardy-littlewood-Sobolev theorem, for p G (1,2], we have L p G 
#" S (IR 3 ) with s = 3(i - |) G [0, §). Then, from Theorem 1.2, the following optimal L p -L 2 
type decay results are obtained. 

Corollary 1.1. (L p -L 2 time-decay estimates) Assume that ||«o||ij"^(R3) < <5o <C 1. If 
further, uq(x) G L p (M 3 ) with some p G (1,2], then the following decay results hold for any 
integer I with < / < N: 

\\D l u(t)\\ H »-i(R3) < Co(l + ty^p- 1 ^. (1.11) 
\\D l q(t)\\ HN+1 - l{m < C (l + t)-t(?-5)-^. (1.12) 

Remark 1.2. All the decay results above are obtained without the smallness of the initial 
perturbation in L p (M. 3 ),p G (1,2] or H~ S (M. 3 ). The results generalize those in [22j[33] for the 
case of three-dimensional space. The similar problem in R n with n = 1, 2 will be investigated 
in future. 

Notations. In this paper, D l with an integer I > stands for the usual any spatial derivatives 
of order I. For 1 < p < oo and an integer m > 0, we use L p and W m ' p denote the usual 
Lebesgue space L p (M n ) and Sobolev spaces VK m,p (R n ) with norms || • \\lp and || • ||vy m 'P, 
respectively, and set H m = W m ' 2 with norm || • \\n m when p = 2. In addition, for s G R, we 
define a pseudo-differential operator A s by 

A s g(x) = [ \H\ s g(0e 27T ^ x< d{;, 

where g denotes the Fourier transform of g. We define the homogeneous Sobolev space H s 
of all g for which ||<?||ij s is finite, where 

Mrs ■= I|a s 5||l2 = Ill£l s 5lli2- 

Throughout this paper, we will use a non-positive index s. For convenience, we will change 
the index to be "— s" with s > 0. C denotes a positive generic (generally large) constant that 
may vary at different places. The integration domain 1R 3 will be always omitted without any 
ambiguity. 
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The rest of this paper is arranged as follows. In the next section, some Sobolev type 
inequalities and some preliminaries are given for later use. Section 3 shows the proof of 
Theorem 1.1 by using Fourier splitting method. In the last section, we obtain the time-decay 
estimates stated in Theorem 1.2. 

2. Preliminaries 

Firstly, we give some Sobolev inequalities which will be used in the next two sections. 
Lemma 2.1. (Gagliardo-Nirenberg's inequality). Let < m,k < I, then we have 

\\D k g\\ LP <C\\D m g\\ 1 £; e \\D l g\\ e L r 

where k satisfies 



!-*=(! 
p n 



q n J \r n 



Lemma 2.2. ([6], Lemma A. 5) Let s > and / > 0, then we have 



\D l g\\ L 2 < CHD'+^H^IMI^, where 9 



l + s + 1 



Lemma 2.3. ([32], Chapter V, Theorem 1) Let < s < n, 1 < p < q < oo, | + ^ = |, then 

||A~ s 5||l<j < C\\g\\ LP . 

Now, when n = 3, we derive an estimate of Lyapunov-type which plays an important 
role in closing the energy estimates at each l-th level in Section 4. 

Proposition 2.1. Let (u,q)(x,t) be a solution to the Cauchy problem (1.1)-(1.2) in M 3 . If 
the assumptions in Proposition 1.1 hold, we have 

d 



dt 



\D l u{t)f m + \\D l+l u{t)f L2 < 0, < I < N - 1. 



(2.1) 



Proof. First, we transform system (1.1) into the following equivalent decoupled system 



u t - Au t -Au = -(1 - A) fj(u) Xj , x£R n , t>0, 
g = —(1 — A)~y x el", t > 0. 



(2.2) 



Multiplying (2.2)i by u, then a direct calculation gives 



(u 2 + |Vii| 2 )t + 2|V-u| 2 -2div 

3 

-^2fj(. u ) u zj\ Vu \ 2 - 



i=i 



i=i 



2 / f'Mndr) 
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Integrating the above equality with respect to x on M. and using Lemma 2.1, we get 
d 



As a result, we find 



t ,Kt)||^ + \\Du(t)\\lz < C\\Vu(t)\\ L ~\\Du(t)\\t 2 . 

A| Kt) ||^ 1 + ||^(t)||| 2 < . (2.3) 

Since fj(u) = 0{u 2 ) when u — > 0, without loss of generality, let fj(u) = u 2 . In terms of 
estimates for the derivatives of the solution u, one can apply D l on (2.2)i, and multiply the 
resulting equality by D l u, then integrating it with respect to x over R 3 , obtaining 

±(\\D l u(t)\\ 2 L2 + \\VD l u(t)\\h) + 2||VlA(i)||i 2 

= -J2 [ D l (uu Xj )D l udx D l V(uu Xj )D l Vudx =: h + h- (2.4) 

j=i J j=i J 

For Ii, by using Holder's inequality and Lemma 2.1, we have 

3 3 



h = -ir[ D l (uu x .)D l udx = -Y\ J2 [ (D k uD l ~ k u x )D l udx 

3=1 JR3 j=1 0<k<l J * 3 

<C^2 \\{D k uD l ~ k+1 u)(t)\\ L ^D l u{t)\\ L6 

0<k<l 

<C^T \\(D k uD l - k+1 u)(t)\\ e\\D l+1 u(t)\\ L2 . (2.5) 



0<k<l 

When k < [|], by using Holder's inequality and Lemma 2.1, we get 
\\{D k uD l ^u){t)\\ L% < C\\D k u{t)\\ VA \\D l ~ k+l u{t)\\ L , 

<C||J3 m «(t)|t*||l^Vt)||^||«(t)||^||D'+ 1 «(0||^^ 

< C(\\D m u(t)\\ L , + \\ U (t)\\L*)\\D l+1 u(t)hz < ^\\D l+1 u(t)\\ L2) 

(2.6) 

where we have used the fact ||u(t)||^jv < CE 2 <C 1 in Proposition 1.1, and m satisfies 

k 1 ,m L ,„ k . ,1 + 1 1. k .„ 

= ( )x( 1 --; ) + ( ) x i • (2-7) 

3 3 y 3 2' v l + V v 3 2 ; l + l v ' 

As a result, since k < [|], we have 

m= 2(/ + l-fc) £[ 2' 1) - 
When > [|] + 1, from Holder's inequality and Lemma 2.1 again, we find 

\\(D k uD l +^ k u)(t)\\ d < C\\D k u(t)\\ L 4D l+1 - k u(t)\\ L3 

<C||ti(t)||^^||^Vt)||^p m u(t)||^||i^Vt)||^^ 

< C(\\D m u(t)\\ L2 + \\u(t)\\ L2 )\\D l+1 u(t)\\ L i < \\\D l+1 u(t)\\ L ^ (2-8) 
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where m is denned by 



l + l-k 1 ,m 1, k ,1 + 1 I, , 

x 7^-T + (-3- " o) x (! 



3 v 3 2' Z + l v 3 2' v Z + l" 



that is, m = ^ G (|, 1] since A; > 
From (2.5), (2.6) and (2.8), we have 



h = -J2 [ D l (uu Xj )D l udx < h\D l+1 u(t)\\ L 2. (2.9) 
The estimate of I2 is absolutely the same to (3.10)-(3.12) except that we replace I by I + 1. 



In fact, we have 



I 2 = - / D , V(uu ay ).D'Vuda: < Jp i+1 «(t)|| L 2. (2.10) 

Thus, from (2.4), (2.9) and (2.10), we get the estimate of Lyaponov-type (2.1). Then, 
we complete the proof of Proposition 2.1. ■ 



3. Decay results with initial perturbation in L 1 (M") 

In this section, we will give optimal decay results by Fourier splitting method introduced 
in \17\ [28] together with energy estimates. Theorem 1.1 will be proved by the following 
lemmas. First, a straightforward application of Fourier splitting method yields an optimal 
L 2 -norm time-decay rates of solutions as follows. 

Proposition 3.1. If the initial data u (x) G L 1 (M n ) n L 2 (M n ) with n > 1, one has 

IK*)llia(R») < C(||^o|Ui(R«) + INHl^r^X 1 +t)~%. 

Proof. First, multiplying (1.7)i by u and summing up them and then integrating over R n , 
we obtain 



1 d 
2cf£ 



/ u 2 dx = / y fj(u)d Xj udx — / u 2 dx + 

jR n jR n . =1 jR n 

= — J u 2 dx+ I u(l — A)~ 1 udx. 



u(l - A)~ x udx 

(3.1) 



By Plancherel theorem, we get 
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We rewrite (3.2) as follows. 



dt 7 K n n + tJ R n J\(,\<y/j V 1 + l^l 2 n + t 

which is multiplied by (n + t) n yields 



d_ 
dt 



(n + t) B / 


|u| 2 d£ 


<{n+ty f ( 









H 2 d£ 



<(n + i)l£(t)|| 2 ) / 

5 ^l<x/f V n + t J 

<C\\u{t)\\l lm {n + t) n -\n + tr^ 

<C\\u \\ 2 Lim (n + t)^-\ (3.3) 
where we have used H'ii(i) ||z,°°(R™) — II^WIIl^r™) an d the fact ||u(t)||x,i < ||^o||l 1 (R' 1 ) i n 

mm®- 

Integrating (3.3) with respect to t, we have 

||«(<)|lia(Rn) < C(||^o||li(r») + ||n || L 2 (R n ) )(n + t)~*. (3.4) 

This proves Proposition 3.1. ■ 

„ n 

In the proof of Lemma 2.1, we used the essential fact that / N fj(u)u x dx in (3.1) is 

equal to 0. However, for the derivatives of u, the Fourier splitting method above can not be 
used directly. The following is the main reason. Similar to prove (3.1), we have 



~/ \D l u\ 2 dx=[ y2D l fj(u) Xj D l udx- [ \D l u\ 2 dx+ [ D l u(l-Ay 1 D l udx. (3.5) 

2 dt 7r™ JR™ jR n jR n 

„ n 

But, the term / \^ D l fj(u) x D l udx in the RHS of (3.5) is not equal to 0. How to control 

this term? Here we use the existence results and the smallness of solutions in H space 
stated in Proposition 1.1. 

Proposition 3.2. Let 1 < n < 4. Suppose ||«o||_ffJV(R«) "C 1. If Uq(x) G L 1 (]R n ) and 
D l u (x) e L 2 (M n ), one has 



\D l u 



ln) < C7(||t*o|Ui (R n ) + ||zAxo|| L 2 (R n))(l + t)-f- z , / = 1, 2, • • • , N. (3.6) 



Proof. Applying D l (1 < Z < iV) on (1.7)i, and multiplying the resulting equality by D'u, 
then integrating it with respect to x over M n , one has 



™/ \D l u\ 2 dx = - [ y2D l fj(u) Xj D l udx- [ \D l u\ 2 dx+ [ D l u(l - A)~ l D l udx . 

(3.7) 
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As mentioned above, to use the Fourier splitting method as Proposition 3.1, we have to get 

n „ 

y^D l fj{u) Xj D l uAx < 8 / \D l u\ 2 dx, (3.8) 

where < #o <S 1- In the following, we only prove (3.8) for 3 < / < N. For I = 1, 2, the proof 
of (3.8) is similar. We omit it here. Without loss of generality, let fj(u) = u 2 , one has 

I D l f 3 (u) Xj D l udx 

= [ uD l u Xj D l udx + I DuD l ~ 1 u Xj D l udx + V / C l k D k uD l+1 ~ k uD l udx 



2<fc<(Z-l) ' 



< 2\\Du(t)\\ LX{Rn) \\D l u(t)\\l HRn) +C Yl 

2<fc<(Z-l) 

<2\\Du(t)\\ L ^ {Rn) \\D l u(t)\\l 2{Rn 



<2\\Du(t)\\ L ^ {Rn) \\D l u(t)\\h^ 
+C\\u(t)\\l^ n) \\D l u(t)\\% {Rn 

<2\\Du(t)\\ L ^ {m \\D l u(t)\\ 2 L2{Rn 

= 2\\D U (t)\\ LX>iRn) \\D l u(t)\\ 2 L2{Rn 



D k uD l+1 - k uD l udx 



+ C||^(t)|| L 4 (R n ) ||Z) / - fc ^.(t)|| i 4 (R n ) || J D i n(t)|| i2(RI , 



+ C 1 1 -a (i ) 1 1 i7 ( ( ^^ 2 } 1 1 ) 1 1 i 1 ^! 1 



+ c|| u (t)||^ e2) ||^n(t)||™y ^ID'uCOH^ 



(3.9) 



where C|. 



Z 

A- 



and 



6i = —± and 6 2 = - 



I 



I 



(3.10) 



By noticing 1 < n < 4 and 2 < k < I — 1, from (3.10), we know < Q\ < 1 and < 6> 2 < 1. 
Then, it follows from (3.9) and the fact ||u(i)||^jv <C 1 that 



/ D'/iW^W^ \\\D l u 
Combining (3.11) and (3.7), we have 



(3.11) 



— / (D l u) 2 dx < - [ (D l u) 2 dx+ [ D l u(l - A^rtudx. (3.12) 

dt J R n J R n J R n 



From Plancherel theorem, one has 



d 
dt 



\C\ 2l u 



d£ < 



\t\ 2l u 



d£ + 



1 



< 



< 



n + 2l 
n + 2l + t 



t 



l£l>> 

2 



L2i n + 21 + t 



+ 

2 



d£ 



\C\ 2l n 



de + 



i 



de+ 



i 



\C\ 2l n 



d£- 



| 5 |<yi±2i 1 + |Ch 



d£ 



n + 2/+t 



l«l<l 



/ n+2l 



\i\ 2l u 



de 

(3.13) 
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We rewrite (3.13) as 



dt 



\Z\ 2l u 



< 



n + 21 + t J M n 
1 t 



\i\ 2l u 



J\Z\<J&™. V 1 + \t\ 2 n + 2l + t 
Consequently, multiplying (3.14) by (n + 21 + t) n+21 , one has 

(n + 2l+t) n+21 [ |e| 
< (n + 2/ + t) n+2Z 



dt 



2/, 



|$|<w^±2T V 1 + l£l 2 n + 2l + t 



\t\ 2l u 



< {n + 2l + t) n+2l \ 



I*I<1 



f 



< ChWllii(Rn)(n + 2/ + t )«+a-i(„ + 2/ + t)" 41 * 



n + 2l + t 

n + 2l 



d£ 
l£| 2/ d£ 



< C||ito|||i (R n)(ra + 2Z + 1) 



n+2l 



(3.14) 



(3.15) 



Integrating (3.15) with respect to i over (0, t), we find 

\\D l u(t)\\ 2 L2(Rn) = |||e| i ^)|li2 (K „ ) <C(||«o|| L i (R n ) + ||Z) i «o|| L 2 (Mn) )(n + 2Z + t)- IL ^. (3.16) 
Then, we get (3.6). ■ 



4. Decay results with initial perturbation in H S (R 3 ) 



This section devotes to the optimal LP-L (R ) decay rates of solutions to (1.1)-(1.2) 
when the initial data is in the negative Sobolev space H~ S {M?) with s € [0, |). 

The following lemma plays a key role in the proof of Theorem 1.2. It shows an energy 
estimate of the solutions in the negative Sobolev space if _s (M 3 ). Namely, we have 
Lemma 4.1. If £q := ||uo||#jv <C 1, for s £ (0, i], we have 

y (|A- s n| 2 + |A- s Vn| 2 ) dx + J |VA" s u| 2 dx 
< C{\\Du(t)f m + \\D 2 u(t)\\ 2 H1 )(\\k- s u(t)\\ L , + ||A- s Vn(t)|| L2 ); 



(4.1) 



and for s G (^, |), we have 



y (|A" s u| 2 + |A" s Vu| 2 ) dx + y |VA- S u| 2 dx 



< c(ik*)h;; 5 \\Du(t)\\h 8 + \\Du(t)\\%S\D 2 u(t)\\}; s )(\\^ s u 



L 2 + ||A~ s Vu 



L ' 2 , 



(4.2) 
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Proof. Applying A s to (2.2) i and multiplying the resulting identity by A s u, and integrat- 
ing over M 3 by parts, we get 

~ J (|A" S ^| 2 + |A" S Vu| 2 )dx + J |VA~ s n| 2 dx 

:= Ji + J 2 . (4.3) 
For Ji, using Holder inequality, Lemma 2.1, Lemma 2.3 and Young's inequality, we have 

3 

h < C||A-^)|| i2 ^||A-7>Mt)||L2 

3 

^CIIA-^llLallVu^ll^llV^III^pViJIllr 



< Cdl^Hi, + \\D 2 U (t)\\liM- S <t)\\L>. (4.4) 

have used the 
Similarly, it holds that 



Here we have used the facts \ + # < 1 and ~ > 6 



J 2 < C(||^«(t)|li a + p^(t)||i 2 )||A-V«(t)|| i2 . (4.5) 



Combining (4.4) and (4.5), we get (4.1). 

Next, we want to prove (4.2). A direct calculation as (4.4) gives 

3 



L 2 + 3 



L 2 



<C7||A-«(*)IU»P«(t)llLa||«(t)ll L | <C\\A- s u(t)\\ L 4Du(t)\\ L 2\\u(t)\\ L 2 2 \\Du 
<C||n(t)||;^|| J Du(t)|||ri|A-Mi)llL- 



(4.6) 

In the same way, one has 



h < C\\Du(t)\\ 8 L 2* \\D 2 u(t)\\h S W A ~ Su 



L 2 ■ 



Thus we complete the proof of Lemma 4.1. ■ 
With Proposition 1.1, Proposition 2.1 and Lemma 4.1 in hand, we are now ready to 
prove Theorem 1.2. 

Proof of Theorem 1.2 for the case of s £ (0, ±]. Define £_ s (t) = \\A~ s u(t)\\ 2 L2 + ||A~ s Vn(t)|| 2 2 . 
Integrating (4.1) with respect to t, we find for s £ (0, |], 

£-.(*) < £-s(0) + C f \\Du{T)f H2 ^/£ljr)dT. (4.7) 
Jo 
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From (2.1), we have the estimates of integrability of ||D"u||^i with respect to t. As a result, 
we have 

£L fl (t) <£- s {0) + C sup y/£- a (T) < C(l + sup V^M)- 

0<r<t 0<r<i 

This yields £- s (t) < Ci with a positive constant Ci, that is 

||A- S u(t)||i 2 + ||A-Mt)||| 2 < ||A- S ^)||i 2 + ||A- s V«(i)||| 2 < d. (4.8) 
This proves (1.8) for s G (0, \]. 

Next, we recall the estimate of Lyaponov-type (2.1) as 

^-\\D l u(t)\\ 2 H1 + \\D l+1 u{t)\\\ 2 < 0, < I < N - 1. 
We may use Lemma 2.2 to have 

\\D l+1 u(t)\\ r.2 > C||A-Vi)ll7J^II£>V^ l|1+ ^ 



From the above inequality and (4.8) we get for each I with < I < N — 1, 

\\D l+1 u(t)\\h > \\\D l+l u{t)\\h + lcc;^\\D l u(t)C + ^ ] 

> C 2 (\\D l+1 u(t)\\h + \\D l u(t)\\l 2 ) 1+ rh, (4.9) 
where C 2 > is a constant. Thus we deduce the following time differential inequality 

^P'uWH^-, +C 2 (|| J D i n(t)||^ JV _0 1+I ^ <0, forZ = 0,l,--- ,JV. 
Integrating this inequality, one gets for some constant C3 > 

\\D l u(t)\\ 2 HN ^ < C 3 (l + t)- (/+s) , for Z = 0, 1, - - - ,iV. (4.10) 
(2.2)2 and (4.10) yield 

\\D l q(t)\\ HN -i-i < C A (l+t)- !± ^ 1 , for I = 0,--- ,N- 1. (4.11) 



Proof of Theorem 1.2 for the case of s E (^, | ). First we can give from what we have proved 
for (1.8)-(1.9) with s = ^ above that the following decay results holds: 

\\D l u(t)\\ 2 HN ^ <C 3 (l + tr l -3, for/ = 0,l,--- ,N. (4.12) 

As a result, and using (4.2) that for s £ (5, |), 

£-.(<) <£-s{0) + C f \\u(T)\\ 8 ^\\Du(T)\\]^'y/n^dr 
Jo 

<Ci + CC 3 / (l + r)-i-fdr- sup 

J0 0<T<t 

<C 5 (1+ sup V'f-.(r)), (4.13) 

0<T<t 
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which yields (1.8) with s £ (\, §). At last, the proof of (1.9) with s E (\, |) can be treated 
as the case of s € (0, g] above. 

Thus, by taking Co = max{Cj}, we complete the proof of Theorem 1.2. 

1< i<5 
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